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Phase-coherent diffusive transport through mesoscopic hybrid superconductor/normal metal tun- 
neling structures is investigated. For a N — s — S two-barrier tunneling system with bulk S and 
N electrodes coupled by a mesoscopic superconducting constriction s, zero- bias conductance and 
non-linear I-V curves are calculated under the assumption that the dwell times of quasiparticles in 
the s region is shorter than inelastic relaxation time. It is shown that the low voltage conductance 
of this system determined by the Andreev reflection processes may exceed the conductance in the 
normal state and its value is very sensitive to the weak pairing interaction of electrons in the s 
region. We show that even weak pairing electron interaction may result in the significant qualitative 
and quantitative change of the conductance temperature dependence with respect to the case of 
structures with the normal mesoscopic region. We calculate the I-V curves and show that they 
depend on the applied voltage in a non-monotonic way, therefore differential conductance becomes 
negative with increasing voltage. Such behavior is due to the voltage dependence of the order pa- 
rameter in the constriction and the phase difference <p between the S and s superconductors. It is 
shown that if the tunneling processes determine the form of the quasiparticle distribution function 
in the s superconductor, the phase <p is stationary at arbitrary voltages. For quasiparticle tunneling 
interferometers in which the mesoscopic superconductor, s, couples the superconductor, S, and the 
normal metal, N, the zero bias conductance, as a function of the phase difference between the S 
electrodes is investigated. It is shown that the amplitude of the conductance oscillations may exceed 
the conductance of this structure in the normal state. 



Phase-coherent transport in mesoscopic superconductor/normal metal (S/N) systems has been an active area 
of research during the last decade fjj|. The interest in the theoretical investigations was stimulated by impressive 
technological advances and by experimental activity in studying various properties of small mesoscopic structures 
P-fTo] . Interesting phenomena in mesoscopic systems are due to the importance of both the phase coherence established 
in the s constriction by the proximity effect and significant departure of quasiparticles from equilibrium. This is 
completely true for two-barrier structures N — s — S with barriers at the interfaces between the N and S electrodes 
connected by a superconducting constriction s of length d. The dimensions of the constriction transverse to the 
current direction are assumed to be small in comparison with the London penetration depth in the S electrode. We 
consider the system with diffusive transport, i.e. we suppose that the mean free path I in the s region is small 
with respect to the constriction length d. Because of non-conservation of the momentum, the interference of normal 
electron wave functions related to reflections from the barriers is not essential. Nevertheless the coherence of different 
(ordinary and Andreev) reflection processes related to the condensate wave function and nonzero order parameter 
A in the superconductor s is very important because the inter-barrier distance d is supposed to be small in comparison 
with y/hD/A, where D = Ivf/3 is the diffusion coefficient. In what follows we assume that transparencies of both 
barriers 2?i.2, (averaged over momentum directions) are small enough to allow the main contribution to the resistance 
of the system to be due to the barrier resistances. Tunneling processes determine the dwell times Tbi,2 = 2?i,2^f/4c£ in 
the s region which are supposed to be shorter than the inelastic relaxation time Ti n in the superconductor s, so that 
the following conditions should be fulfilled 



where Tdif — fo/(D/d 2 ) is the diffusion time of quasiparticles through the length d . It is clear that the quantum 
nature of the tunneling processes becomes more pronounced if the tunneling rates h/rbifi are comparable with the 
characteristic scale of the quasiparticle energy, 
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fi/ni,2 ~ A, 



(2) 
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because under this condition the classical notion of a quasiparticle whose dwell time should be longer than its energy, 
loses its sense. Nevertheless the Green's function approach enables one to obtain the quantum kinetic equations as 
given in which are valid beyond the classical limits, i.e. when the quasiparticle energy is not large compared to the 
tunneling rate. Note that under the conditions (|l|) the proximity effect, i.e. the influence of S and N electrodes on the 
condensate wave function and on the order parameter in the s region, is strong. We also note that unusual features of 
transport properties are due to the significant role of Andreev reflection processeses in the considered system. These 
processeses occur in the presence of two potential barriers (at x = and x — d) and the superconducting order 
parameter which has a two-step form: A(x) = A8(x)8(d — x) + As 8(x — d) exp(iip), where ip is the phase difference 
between the superconductors arising at non-zero voltage, V and 0(x), is the Heavyside function. As a result, the 
energy dependent transmission coefficient 2? e (A, ip) of quasiparticles with energy e < A<j, that determines the current 
at low temperatures, appears to be a function strongly dependent on V through the voltage dependence of A and 
ip. All these circumstances result in nontrivial features of the quasiparticle phase coherent transport through the 
N — s — S system which will be investigated in this paper, ft should be noted that some of these phenomena have 
been studied in fl5|| . We therefore investigate the transport phenomena in more detail with emphasis upon the case of a 
weak pairing electron interaction in the s region, i.e. when the critical temperature, T c q of the superconductor, s in the 
absence of the pair-breaking and proximity effect (t^i 2 = 1S small in comparison with the critical temperature of 
the S electrode, T c s- It will be shown that in spite of a small ratio t c = T c q/T c s the properties of the N — s~ S system 
may radically differ from properties of a two barrier N — N — S structure well studied in the limit t c = fl4 21|. We 
also study the zero bias conductance as a function of the phase difference between the S electrodes in a quasiparticle 
tunneling interferometer with two superconducting electrodes coupled by a mesoscopic superconductor s. It is shown 
that the amplitude of the conductance oscillation may exceed the conductance of this structure in the normal state. 



B. The N-s-S system 

We consider the N — s — S system shown in Fig. la. As in Refs. Jf4|-p2[ we use the approach based on the equations 
for the quasiclassical Green's function G — G(r,p F ; e) which is the 4x4 supermatrix(see Ref. |TT 

G= {6 G A 

consisting of the retarded G R , advanced G A , and Keldysh G K , Green's functions which are 2x2 matrices in Nambu 
space. Note that we suppose that a stationary solution realizes and the Green's functions do not depend on time. 
This non-obvious assumption is justified by the final result. The matrix G K is related to the matrix distribution 
function f = f 1 +f& z 

G K = G R f - fG A (3) 

The matrices G R ' A have the following form 

G' 1 = g^ z + />, /> = f»ia v exp^x) 
where \ is f ne phase of the order parameter /1 = R{A). The current in the system is given by the following relation 

/ = — -j-Tr<7 z / de{G R d x G K + G K d x G A ) (4) 



where A = w y w z is the cross-section area of the s region. The transverse dimensions w VtZ should be small compared to 
the London penetration depth. Therefore we need to solve a one-dimensional equation in the s region (0 < x < d, x- 
axis coincides with the direction of the current), where in the considered diffusive case the matrix G = G(x, e) averaged 
over the momentum direction obeys the equation (see Ref. [pT[) 

Dd x (Gd x G) + i[ea z + A, G] = 0. (5) 

a z = \a z is the Pauli supermatrix, and the order parameter supermatrix in the film is A = lA where 



A = 



A 

-A* 
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The order parameter is given by the self-consistency relation which in the framework of the weak coupling theory has 
the form 



A = A 



de(f R f-fn 



(6) 



where the constant A determines the critical temperature T c q of the superconductor s in the absence of pair-breaking 
factors and the proximity effect, 

T c0 = 1.14w D exp(-l/A) 
The matrix G obeys the normalization condition 



G 2 = l 



(7) 



In order to solve Eq.(||) we need to take into account the boundary conditions |T^] that in the diffusive case reduce 
to pi (see also @). 



D(Gd x G)(+0) = e bl d[G(+0),G N ], D(Gd x G)(d-0) = e b2 d[G s ,G(d-0)] 



(8) 



where ey = pD/2dRbjn, Rbi,2o is the interface resistance per unit area at the N/s (x = 0) and s/S (x = d) interfaces, 
Gs,n equilibrium Green's functions in the electrodes, p is the normal state specific resistivity of the superconductor s . 
Note that the energies, tbj , are connected with the characteristic dwell times: Tbj = hjtbj- In terms of the Thouless 
energy D/d 2 = Er^the conditions (jj]) may be written 



h/T in < €bj < E T h 



(9) 



Suppose that the length d of the s region is small enough i.e. d <C \ffiD / As- Then the solution of Eq.(||) is readily 
found(see Appendix I). The retarded and advanced Green's functions are given by 



G»(e)=g»(e)a z +F(e) 



eK(e)a z + [Aia y +ie b 2fg(e) 



(10) 



where 



C(e) = [(^{e)f - A 2 - 2Aie b2 /£cos^ + (e &2 /£) 2 ] 1/2 



(e) = e + ie b 2gg' A {e) ± ie bl 



The Keldysh function is given by Eq.(A8) and it is convenient to separate the anomalous part G^ to give 

G K = G R n - nG A + G*f 



We have for G* 



with the anomalous self-energy 



where n(e) = tanh(e/2T), 



Gf = {Ef-G«EfG A ) 



1 



(C R + ( A ) 

= 2ie 61 {n_(e) + a z [n+{e) - n(e)]} 
n±{e) = [n{e + eV) ± n(e - eV)]/2 



(11) 



(12) 



(13) 



Thus the anomalous part G^ is determined by E^ which contains only the self-energy depending on the N-electrodc 
Green's function G K . Using Eq.(0) enables one to find the non-equilibrium part of the matrix distribution function 



5f = (fo-n)+fa z 



(14) 



which determines the anomalous part of G K : 
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Gf = G R Sf - SfG A (15) 
From Eq. ( |ll| ) we find for the non equilibrium parts of the distribution functions 

f= Av{C R + QA) ™«{l-G A G«) (16) 

Sfo = m^T^) Tt ^ {(7z ~ dA&zdR) (17) 

where v — Heg R is the density of states, Sfo — fo — n,. Using Eq.(|l2|) for the non-equilibrium part of the 
distribution functions may be written in the following form 

Sfo = a+(n+ — n) + bn~ (18) 
/ = a_n_ — b(n + — n) 



where 



e bl M± Ae fcl e fc2 Re/f . 

a± — - — - — —5 , b — — — S- sm w 

2vlmC, R i/ImC fl |£«| 2 



M± = 1 - g R g A ± [A 2 - 2e fc2 A Im /# cos p + (e b2 \f R \) 2 } ' 



IC*| 2 



We took into account that ( A — — (C R )* , g A — ~(g R )* . From the self-consistency relation (ph we obtain the following 
system of equations for A and ip, 

A A = €b2 {a cos ip — (3i sin ip) , (19) 
/3 A = £{,2 (a sin 93 + /3i cos^), (20) 



where 



A = \n(T/T c0 ) — J de(f (e)Re 



1 n(e) 



C«(e) e 



(i 



&/ (e)Im^| 



A = r def(e) Im - ^M^^ - fc = °> L 
Jo C K (e) 

Note that Eq. ( pp| ) is the consequence of the current conservation law in the x-direction. Introducing the normalized 
order parameter 5 = A/e^ one can reduce Eqs. ( |l9| ) and (|2(]) to the equivalent ones 



v oA+^i+i(a^o-A^i) , . 

exp(i^) = 7 ===- T === — (22) 

From the self-consistency equations (jl^) and (^0|) or (|l]) and (p2|), it follows that the stationary solution for the 
order parameter exists at arbitrary V and transition to the ac Josephson effect (time dependent phase difference 
ip) does not occur with increasing voltage. In other words the critical current of the S/s tunnel junction is absent 
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in the considered mesoscopic system. Such a situation differs radically from that in a single S/S tunnel junction 
composed of two bulk superconductors. If at least one of the two superconductors has mesoscopic dimensions, it is 
important how it is connected with the conductors and non-equilibrium states arising in the presence of the current 
play a significant role in this case. In our system one of the important aspects of the non equilibrium state is 
the quasiparticle charge-imbalance determined by the distribution function /(e) , and, as a consequence, the gauge- 
invariant potential fi = $ + (h/2e)dtX m the s region, where 4> is electrical potential and \ is the order parameter 
phase. Under the assumption ([!]) the solution for the phase difference between the superconductors is stationary for 
arbitrary voltages. Therefore we can set x = in the s region so that ^ — \ / °° def(e)i/(e) coincides with the voltage 
between the superconductors $ ^ —(h/2e)dt<p = 0. In other words the Josephson relation between the frequency 
(equal to zero) and the voltage drop across the superconducting tunnel junction is violated in the structure under 
consideration. The current may be calculated at any point x and, for example, at the N/s interface we obtain 



2eR 



N J- 



de{F_(e)n_(e) + F+(e)[n+(e) - n(e)]} , (23) 



where 



F_(e) = (1 + r)u{e)[l - a_(c)], F + (e) = v(e)b{e){l + r) , 

r = Rb2/Rbl = £bl/e&2 ■ 



If the S electrode is a conventional BCS superconductor, 



g§(e) = fi(e)e/A s = e/y/(e + i0)* - A% . 
Then for |e| < As, b(e) = 0, F + (e) — and at \eV\ < A at zero temperature the current reads 

/ = -L fdeF^e). (24) 
eHN Jo 

Note that the function F_ (e) = F- (e; V) depends on voltage through the voltage dependence of A and ip. It represents 
the transmission coefficient of the system which determines the efficiency of Andreev reflection processes. Taking into 
account that 6(e) = and assuming |eV| < As , we find for the non-equilibrium part of the distribution functions 

/(e)=a_(e)sgn(eV)0(|eVl-M) > (25) 
Sf (e) = -a + (e)sgn(e)6(\eV\ - |e|) . (26) 

Consider the case of small critical temperatures of the superconductor s, T c q/T c s 1, and also assume that the 
following condition is fulfilled 

A, e b i, e b2 ,eV < A s . (27) 
In this case Eq.([2l|) and ( p2| ) for A and cp can be simplified and presented in the form (see Appendix II). 




A" + /38' (28) 



cos ip + 1 sin ip 



From Eq. 



and (A10) we find the current 
(e&i + e h2 )/ei?Ar 



ak- fil +i/3 (a + A) 



du 



\C(u,n v )\' 



(29) 



(30) 



The I-V curves obtained by numerical calculations results on the basis of Eqs.(Eq), ( p9| ) and (30) are presented in Fig. 2. 
One can see that (as a consequence of the order parameter suppression in the s region) the differential conductance 
becomes negative with growing voltage. 
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In general the solution of Eqs.(E8h, (|29|) and the current can only be determined numerically because the formulas 
are rather complicated. Nevertheless the zero-bias conductance go = G(0)/Gn can be found from Eq. (|30|) , where 
G(V) = dl/dV . It is given by 



9o(5, r) 



{l + r)r(l+S) 2 



[r 2 + (1 + S) 2 f/ 2 

where according to Eqs.(28) and ( |l9|) S = A/e^ is defined by the equation 



So t c 



(31) 



(32) 



It follows from Eq.(B2|) that (under the condition (|27|)) the proximity effect is strong, i.e. A 3> Ao. Moreover 



A 

— ► oo at Ao 

Ao 







i.e. due to the proximity effect, anomalously big enhancements of the order parameter occur at very weak pairing 
electron interaction in the s region. Assuming S -C 1, one can obtain from (p2[) that 



f 



In 



4A f 



e b2 (r+Vr J +l) 



A So ln (r+V^+T) 

This expression is valid for very small So which satisfies the condition 

4A S 



(33) 



, (r + Vr 2 + 1) 
In ^ '- > ln 

So 



e b2 {r + Vr 2 + 1) 



that is fulfilled provided So <C (ebi + £f,2) 2 /A| <C 1. In particular if e b \ + eb2 ~ 10~ 2 As the requirement So <C 
I0~ 4 means that (B3|) is valid provided T c q is anomalously small, T c o <C 10~ 6 T c s, then A ^> I0 4 Ao. It can be seen 
from condition (|27) that one can ignore the presence of the order parameter in the s region (i5 <C 1) only if the pairing 
interaction in the s region is very weak, i.e. T c o <C e^efci + £b2) 2 /A| . The zero-bias conductance as a function 
of t c is shown in Fig. 3 for different parameters r = Rb2/Rbi (and tb2 = O.O5A5). One can see that the normalized 
conductance may be both smaller and bigger than unity. In particular from (Bl]) we find that for r > l/\/2 the 
maximum value of the conductance, g corresponds to S — S rn where (I + 5 m ) = \2r and Eq.(|T]) equals 

(1 + r) . 



3V3 



From Eq.(|32j) we find that the maximum conductance is realized for the case when the critical temperature T c o — X^,™, 
where 



T. 



AT, 



cS 



C6l(l + V3) 



AT, 



cS 



l/(l-l/v/2r) 



(34) 



Eq.(|34|) is applicable for r satisfying the condition T c ™ <C T c s; in particular it is true for r > 3V3/2 — 1 which 
corresponds to 50 max > 1- At low temperatures T <C A5 for zero-bias conductance we find from (E3I) 



g(t) = 2(l + r)fi(f) 



du 



Re(2tu + ir)/C(2tu,n(t)) 



la cosh 2 u {2tu) 2 + r 2 + fi(i) + \Q{2tu, ft(i))| 2 
where t — T/e b2 , Q(t) = (1 + S(t)) 2 and £(u, fi) = [(it + ir) 2 — SI] 1 / 2 , the function S(t) is defined by the equation 



(35) 



(36) 



with 
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4A & 

a(Q,t) = In / du -Re- 



(Vti + r 2 +r) e b2 Jo expu+1 £(tu,Q) 



A(n,t)=ln h / j - ln — r 

oo Jo cosh u y/fl + r 2 + r 



The results of numerical calculations on the basis of Eqs.(|35j), Q3q ) are presented in Figs. 4 and 5. We see that 
the conductance may be a non-monotonic function of temperature that radically differs from the corresponding 
dependencies occurring in the case of normal mesoscopic region with T c q — shown by dashed lines in Figs. 4 and 
5. Thus a weak pairing electron interaction results in significant qualitative (for r > 1) and quantitative changes of 
the conductance dependence with respect to the case of structure within the normal mesoscopic region. Fig. 6 shows 
that if r < 1, the conductance may be non-monotonic function of temperature even at T c q — 0. We see that the 
pairing electron interaction results in a shift of the position of the conductance maximum to higher temperatures 
together with an increase in width of the maximum. The latter is due to the slow decrease of the order parameter 
with increasing temperature. 

Consider the case when the resistance of the barrier at the N/s interface is small enough (r » 1). To be more 
exact we suppose that 

e 6 i > A , e b2 . (37) 

i.e. r > A In this case the energy gap is absent in the superconductor s due to a strong pair-breaking effect of the 
normal electrode. If the condensate Green's function is small, all the expressions are significantly simplified and from 
@ at T, eV < A s we find for the current 

/ = (ebl J" eh2] (S 2 + 25 cos tp + 1) Im tt(r + teV/2nT) , (38) 
eitjv 



where ^(2) is the digamma-function, T = 1/2 + e&i/27rT, 5 and costp are given by Eqs.(|28j), ( |29| ) with 

A = hx(T/T c0 ) + Re^(r + ieV/2nT) — *(l/2), /3 = Im *(r + ieV/2i:T), (39) 



a = A s / de- 



(e 2 + e 2 bl )^A 



2 _ e 2 



At zero temperature we obtain from Eqs.(|38|), ( p9| ) 

I 



(e&i + £b2)/ei?Ar 



(5 Z + 2(5cos<y9 + 1) arctan(v/r) , (40) 



where Eqs.(39) reduce to 



. . 2^/v 2 + r' 2 25 D , . , 

A = ln , a = ln , f3 = arctan(w/r) . 

The TV curves computed for this case are shown in Fig. 7. At small voltages we have IRn = g^V , where the 
normalized conductance is given by the expression 



.9o 



In 2 j- 

00 



At large voltages, e b i "C eV <C A5, the normalized current has the form 

1 \n 2 {2A s /eV) + ln 2 (2 e y/A ) + 2y / [ln 2 (2A s /eV r ) + l][ln 2 (2el//A ) + 1] 

(cm + e b2 )/eR N ~ ln 2 (2ey/A ) + 1 

At (2eV) 2 > AsAq this function slowly decreases with increasing voltage. 
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C. Quasiparticle interferometer 



Consider a quasiparticle interferometer composed of three tunnel junctions (see Fig. lb) in which the phase difference 
f between two different S/s interfaces is set by an external magnetic field. A similar system in which S and N 
electrodes were in contact with a normal metal was considered in fl^-|is||. Suppose that two barriers at the S/s 
interfaces are symmetrical with resistances equal to Rb2 and the resistance of the barrier at the N / s interface equals 
Rbi ■ We again assume that the resistance of the system is determined by the barriers and in the normal state is given 
by the expression Rn = Rb2/2 + Rbl- Assuming that the width of t he s uperconductor s is small, W <C y^D/A one 
can neglect the spatial variation of the Green's function. Then Eq.(Al) is valid with 

£ = ie b2 Gs+ + itb2Gs- + efciGjv , (41) 

where the Green's functions Gs± correspond to the phases ±93/2, e^j = pDwj /2dWRbjn , w\ = W being the width 
of the s region and W2 is the width of the S/s interfaces. As in previous cases we find 

^ (e )=^ + /M( e) = !!^^^ l (42) 

where 

^{e) = {{e»{e)f-[^{e)f} 1 ' 2 , 



c 



R - A (e) = e + 2ie b2 gf A {e) ± ie bl , A"(e) = A + 2ie &2 /£(e) cos^/2) 



From the self-consistency relation ^ at zero voltage between the S and N electrodes the following system of equations 
for A and if can be found 

roo f R (e) 

AA = a = -2e b2 cos(y>/2) / den(e) Im ^IrrT • ( 43 ) 

Jo C (ej 



where 



A = ln(T/T c0 ) - J™ de (rc JL- - l\ n{e) 
At T = assuming as before that Aq, 0,1,2 -C Ag , we obtain 



A2+ £ 2 + £w 

A = ln^ -r , (44) 

^0 

4A S 

a = 2e b2 003(95/2) In 5 , (45) 

where A^ = A + 2e b 2 cos(<^/2). It is convenient to introduce the function <L : A/2eb2 — 5 v cos(ip/2), then from 
Eqs.(E3) and (MS) the following equation for S v can be found 



v /(^ + l) 2 cos 2 (^/2)+r 2 +r 4 



(^ + l)ln^ = ln T' ( 46 ) 

00 t c 



where Sq — A /2ef, 2 - After calculations similar to those carried out in Refs. |16|-{19f we obtain for the zero-bias 
conductance of a symmetrical quasiparticle interferometer at zero temperature 

G(0,<p) _ (l + r)r(l+^) 2 cos 2 (<^/2) 



G N [r 2 + cos%/2)(l + <5 v ) 2 ] 3 / 2 



(47) 



Note that at tp = Eq. ( f47| ) is identical to Eq. ( |3l| ) . Thus the amplitude of the conductance oscillations may exceed 
Gn if the ratio r is large enough. For r 2 ^> (5 + l) 2 (when the energy gap is absent in the s region), Eq.MTn yields 



«^-(*). (48) 
<~tn rln (2r/d ) 

Since 4/t c ^> 2t/8q 3> 1 (Ag 3> e&i) the amplitude of the conductance oscillations appears to be much larger than in 
the case of interferometer with a normal mesoscopic region ]l6| |l9| 
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D. Conclusion 



In conclusion, we have studied phase-coherent diffusive transport through different tunnel structures with S and 
N electrodes coupled by a mesoscopic superconductor s. Our study has centered upon the case of a weak pairing 
electron interaction in the s region which defines a critical temperature T c0 <C T c s- If the dwell time in the s region, 
determined by the tunneling processes rt, is small or comparable with h/Ao, the proximity effect is strong, i.e. the 
order parameter A 3> Ao ~ T c q . As a consequence, the subgap conductance of an N — s — S tunneling structure and of 
a quasiparticle tunneling interferometer, depends strongly upon the pairing electron interaction in the s region when 
Too <C T c s ■ Depending upon the ratio of the barrier resistances, the value of the subgap conductance, determined by 
Andreev reflection processes may be both larger and smaller than the conductance of these structures in the normal 
state. We have shown that even weak pairing electron interaction may result in the significant qualitative (in the case 
r > 1) and quantitative change of the conductance temperature dependence with respect to the case of structures 
with the normal mesoscopic region. The subgap current non-monotonously depends upon the voltage, due to the 
suppression of the order parameter in the mesoscopic superconductor. 
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I. APPENDIX 



Integrating Eq.(g) over x and taking into account the boundary conditions, we obtain the following equation for 
G (in the following G denotes the function averaged over the length d ) 

[E,G]=0, (Al) 

where 

E = ea z + A + E, 
S = ie bl G N + ie b2 G s (A2) 



When writing equations (||) and (Al) we disregarded inelastic collisions due to condition ([!]). Let the potential of the 



superconducting electrode be zero and the potential of the normal electrode be equal to V, so that 

G R ' A = ±a z , G% = (1 + a z )n{e + eV) + (1 - a z )n{e - eV), (A3) 

G R S A = 9f A o z + f*' A , Gf = (G§ - G A s )n{e). (A4) 



When finding the solution of ( |A1| ) it is convenient to let the phase of the order parameter in the s layer equal zero 
and the phase of the superconducting electrode S be equal to the phase difference tp which arises in the presence of 
the current. Then from dAl| ) and ([?]) we find the expressions for the retarded and advanced Green's given in Eq. ([To|) . 
The equation for G K has the form 

E R G K - G K E A = G R E K - E K G A (A5) 
where E K = Ef + E^. It is useful to take into account that 

E R > A = ^AqR,a (A6) 

Then using Eq.(Q), we have 

G R G K + G K G A = 6 (A7) 

Therefore from (|a|) it follows that E R G K - G K E A = (( R + <^ A )G R G K and from (|a|) we find for the Keldysh 
functions 

G K = (E K - G R E K G A ) -n^ x - ) (A8) 
II. APPENDIX 

Here we present simplified formulas for f2 ,a and by taking into account the following identity which may be 
readily proved for small energies e< Aj : 

(2^1mC a )(u,n y )=r (u ' + r3+ y +r, 
\([u,%)\ 

Using the notations u = e/efc2, (,{u, 0) = [{u + ir) 2 — fi] 1 / 2 , fl^ = S 2 + 25 cos <p + l one can obtain the expressions for 
a± from M). 



a+ = 1, (A9) 
w 2 +r 2 + fVHC(u,^)r 



°- = 1 - - — 7777 ^ .,2 - (A10) 
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At zero temperature and eV <C A$, one has 



A = A - / deSf (e)Re-^, (All) 



where 



A = In 



r 2 + r 



<5o 



5 = A /e 6 2, A = 1.76T c0 , 



Ao being the gap of the superconductor s at T = in the absence of pair-breaking factors and the proximity effect 
(ey = 0). Introducing the normalized voltage v — Ve/eb2, from Eq.( All ) we find expressions for a and (3 , (see 
Eqs.p| and (ET 



A = In 



|C(^,^ y ) + ^ + ir| 
<5o 



(A12) 



a = In 



4<5 



2 Jo f^) 
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s 



N 

(a) 



FIG. 1. The N/s/S system under consideration, (a). The S/s and the N/s 
interfaces have Rb2 and Rbi barrier resistances respectivly. 
(b). The schematical representation of the Andrccv interferometer. 
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FIG. 2. I-V curves at zero temperature for r less than 1, 662/Ag = 0.05. 
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0.00 0.02 0.04 0.06 0.08 0.10 



FIG. 3. Dependence of the normalised conductance, g on t c . 




O.O 1.0 2.0 3.0 4.0 5.0 6.0 7.0 



T/e B2 

FIG. 4. Temperature dependence of the zero-bias conductance of the N — s — S structure for 
t c = 0.0 and 0.05, e o2 /A s = 0.05 : r = l,r = 2,r = 3. 
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FIG. 5. Temperature dependence of the zero-bias conductance of the N — s — S structure for 
r = 1, e o2 /A s = 0.05 : t c = 0.0 , 0.01 , 0.02 , 0.04. 
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FIG. 7. I-V curves at zero temperature for r> 1, £62 /As = 0.05. 
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